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Abstract 

A novel DNA molecule stretching technique is developed and tested herein. Through a heated converging- 
diverging nnicrochannel, thermal convection and thermophoresis induced by regional heating are shown to 
significantly elongate single DNA molecules; they are visualized via a confocal laser scanning microscopy. In 
addition, electrophoretic stretching is also implemented to examine the hybrid effect on the conformation and 
dynamics of single DNA molecules. The physical properties of the DNA molecules are secured via experimental 
measurements. 
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Background 

The past two decades has witnessed a tremendous growth 
in knowledge regarding the mechanical properties of 
DNA and its polymeric behavior. In addition, develop- 
ments in molecular biology and micro- or nanotechnology 
have increased the interest of scientists and engineers in 
the mechanical manipulation of single DNA molecules. In 
fact, engineering DNA stretching would be a key step in 
the development of the next generation of biological 
microfluidic devices [1]. 

The ability to directly manipulate and visualize single 
DNA molecules has led to a number of advances in our 
current understanding of the physical and biological prop- 
erties of DNA. Two general approaches to DNA stretching 
are in common use: DNA is stretched in a solution as it 
flows through a microchannel or it is stretched on a solid 
surface. Both approaches have their own advantages/disad- 
vantages which depend on the particular application. For 
the former, with fluorescently labeled DNA molecules, it is 
possible to visualize the change in the conformation of a 
single DNA molecule under an optical microscope [2,3]. 

Recently, Ichikawa et al. [4] have presented a novel 
DNA extension technique via laser heating. They proved 
that the new stretching technique was promising and 
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could work in selected applications. Thermophoresis has 
also been found to play an important role in DNA mol- 
ecule stretching. 

The thermal convection induced in this study was 
similar to the convection that is inferred for the well- 
known Earth's mantle convection/or Bernard cell con- 
vection. Such convection produced the horizontal flow 
which caused the movement of the solution. Following 
[4], the governing equations of thermal convection in 
the study are the conservation equations of mass, mo- 
mentum, and energy with the major dimensionless par- 
ameter of the Rayleigh number, indicating the vigor of 
convection and nondimensionalized heat flux. 

In this study, a more powerful modified method was 
used to stretch a single DNA molecule in an electro- 
osmotic flow (EOF) through a heated converging-diverging 
microchannel, using both confocal laser scanning micros- 
copy (CLSM) and micro-particle image velocimetry ((iPIV) 
measuring/visualization techniques. Figure 1 shows sche- 
matically the gradual contraction (8:1) /gradual expansion 
(1:8) flow cell system used in this study. Our main focus 
was to examine the contribution of stretching due to ther- 
mal convection, thermophoresis, electrophoresis, or a 
combination thereof in order to gain further insights into 
the flow behavior of the DNA stretching mechanism and 
the physical/mechanical properties of single DNA mole- 
cules, as well as related phenomena. 
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Figure 1 MicroChannel geometry and observed sections. 



Methods 

PDMS flow cell fabrication 

For this study, we used a 400 x 50 [im and 50 x 400 [im 
converging-diverging test section with a heating foil, 
which is a silicon-based heater with a size of 20 x 5 x 2 
mm, with a total electrical resistance of 20 H, connected 
to a direct current (DC) power supply (N6731B DC 
power supply module) embedded underneath the back- 
side of the floor of the channel The size and dimensions 
of the heating foil were chosen and designed so that the 
temperature distribution on the xz plane (at y = 0) of 
the test section remained uniform upon heating. The 
microfabrication process followed that of [3], except for 
slight modifications in the channel size and converging- 
diverging ratio. The relevant geometric size and dimen- 
sions are listed in Table 1. After completing (8:1:8) the 
fabrication, the test channels were rinsed in acetone 
and ethanol and dried with an argon stream. The 
present study used untreated/treated polydimethylsi- 
loxane (PDMS) channel to measure electrophoresis 
(DNA molecules) velocity and total velocity of EOF, 
respectively. 

For the tested channels, precise information on the 
channel dimensions was extremely important in order to 
make an accurate evaluation. The depth, width, and length 
were measured optically within an accuracy of ±0.2%. To 



understand the surface condition of the present device, 
the roughness of the channel was measured along its cen- 
ter with a surface profilometer. 

Each of the reservoirs (up/downstream plenum) had a 
volume of 0.15 ml. The channel had a total length of 30 
mm, with a length of 800 \im for the test section. The 
detailed values of the test cells are listed in Table 1. 

CLSM/mPIV and mLIF setup 

The CLSM measurement setup, as shown in Figure 2, is 
combined with a laser light source (Ar-ion laser 488 
nm/ HeNe laser 532 nm) and scanning system in order 
to generate the entire field. The flow cell was mounted 
onto an epifluorescent microscope (IX71/FV300, Olym- 
pus, Tokyo, Japan) equipped with a x40 magnification, 
NA 0.85 air immersion objective lens, following that de- 
scribed by [3]. The EOF was driven by a high-voltage 
power supply (PS 350, Stanford Research System, Sunny- 
vale, CA, USA) to drive the flow, with a slight modifica- 
tion for the flow cell and the flow circulation loop. For 
that reason, all the details have not been repeated here. 
The experimental scheme used to implement the (iPIV 
measurement is shown in Figure 3. The use of the (iPIV 
technique is very attractive in microfluidics because it 
helps to determine the detailed flow phenomena of 
microsystems by utilizing flow-tracing particles to map 
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Table 1 Relevant parameters 



Parameters 



Value 



Channel total length, U 
Channel test section length, Ls 
Channel contraction length, 
Channel main width, 
Channel contraction width, Wc 
Channel depth, H 
Channel hydraulic diameter, Dh 
Channel contraction ratio 
Channel expansion ratio 
Electric field (kV/m), 
DNA concentration, jjg/ml 
Working fluid 
Viscosity (cP), \i 
Reynolds number. Re 

X-DNA contour length (pm) (labeled with YOYO-1) 

Radius of X-DNA gyration (pm) 

Temperature ( C), T 

Relaxation time (s), ir (Rouse model) 

Relaxation time (s), ig (Experiment) 

Deborah number 

Velocity vector distribution 



25 
0.0456 



30 mm 
0.66 mm 
0.2 mm 
0.4 mm 
0.05 mm 
0.1 mm 
66.67-160 Mm 

8:1 

1:8 
5, 7.5, 10 
0.065 
IxTBE 

1 cP 
0.032 ~ 0.064 

21 

0.7 



35 
0.0441 



45 
0.0427 



0.6 
1 .2 ~ 2.3 



55 
0.0414 



W ^ 

lOOnm/s 



the flow in the microchannels. In this study, the stained 
DNA molecules could also be used as seeding. 

The setup shown in Figure 3 was based on two pulsed 
Nd:YAG lasers (New Wave SoloII, New Wave Research, 
Fremont, CA, USA; 30 mj, double cavity) firing on the 
second harmonic SoloII (green, 532 nm). The laser provided 
a laser beam with a measured area. The light was positioned 
so as to illuminate the entire inlet, outlet, and midsection of 
the channel. The laser pulse duration was 4 to 80 ms, based 
on the velocity magnitude. The test system was mounted on 
a movable xz stage on an inverted epifluorescence micro- 
scope (DMILM, Leica, Solms, Germany) with xlO magnifi- 
cation, 0.25-numerical aperture panchromatic objective, and 



a field view of 800 x 600 \im. The measurement plane 
(i.e., the object plane) was precisely positioned relative 
to the test section by vertically moving the objective 
lens in the y direction and by horizontally moving the 
table in the x and z directions. The concentration of 
stained DNA molecules based on the interrogation vol- 
ume was 8 X 10^ particles/ml. 

The images were recorded using a Dantec 80C77 
HiSense PIV (Dantec Dynamics, Ulm, Germany) 1,344 x 
1,024 X 12 bit interline transfer camera. Five images 
were taken for each flow field, with a spatial resolution 
of 32 X 32 pixels. The interrogation cell overlay was 50%. 
Background-noise influence was removed by subtracting 
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Figure 3 Schematic of the jjPlV/iaser-induced fluorescence (jjLlF) system velocity and concentration measurements. 
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the background intensity from the captured images. In 
addition, an ensemble averaging 20 images consecutively 
captured for 4 s was used to obtain the velocity measure- 
ments. The calculated measurement depth of the present 
(iPIV was 50 [im. Each measurement was repeated at least 
three times under specified conditions. The measurements 
were conducted in the middle region at both the inlet and 
exit regions of the microchannel. The flow was found to 
have reached full hydrodynamic development at the mid- 
dle region of the microchannel. 

Visualization of the local buffer solution temperature 
was achieved with the same apparatus used for flow 
visualization and measurements (see Figure 3). However, 
instead of using stained DNA molecules, the channel 
was filled with a solution of rhodamine B, a fluorescent 
dye which shows a temperature-sensitive quantum yield 
in the range of 0°C to 100°C [5,6]. Experiments were 
conducted with a fluorescence microscope equipped 
with a long-working distance xlO objective lens. The im- 
ages were recorded with the same equipment used for 
the (iPIV measurements. From the captured images, the 
detailed temperature distribution could be extracted. 
Following [5], the intensity values of the captured images 
were converted to temperature using intensity-versus 
-temperature calibration; calibration of the intensity of 
temperature was made for each solution. 

Flow system 

In the electro-osmotically driven flows, a 30-mm-long 
converging (8:1) -diverging (1:8) microchannel with a 
cross section of 100 x 400 [im and two reservoirs (up/ 
downstream plenum) was used to supply a buffer of 
stained DNA molecules for the channel. Before use, the 
microchannel and entire flow loop were rinsed with DI 
water for at least 1 h to remove any contaminants. The 
transparent nature of the microchannel surfaces allowed 
visual examination of the channels to ensure that no 
bubbles were left. The buffer solution used was Ix Tris- 
borate with ethylenediaminetetraacetic acid (EDTA) 
(TBE) with pH 8.3. A schematic diagram showing the 
flow cell and the auxiliary system is given in Figure 3. 

During each measurement, the microchannel was 
connected to small reservoirs. Current data were recorded 
from the power source by a personal computer-based data 
acquisition system. (iPIV measurements were tal<en 
through a viewing window at midplane {y = 0) between the 
two cylindrical reservoirs with a diameter of 5 mm. The 
potential was applied via platinum electrodes immersed in 
the two 0.15-ml open reservoirs. The distance between the 
two reservoirs was 30 mm. When electric field was >10 
kV/m, the EOF velocity of the solution will increase, and 
the mobility would be dependent on the electric strength 
[6,7]. In order to avoid joule heating, electric field strengths 
of 5, 7.5, and 10 kV/m were thus applied. 



The (iPIV measurement system included visualization 
and the capture of images, the calculation of two- 
dimensional velocity vectors, and post-processing for 
data analysis. The vector field of the flow velocity within 
the measurement plane of the light sheet was deter- 
mined by measuring the displacement of the tracer par- 
ticles and the time durations of two laser pulses. A PIV 
2100 processor (Dantec Inc.) was optimized to process 
the (iPIV images into a raw vector map in real time and 
to transfer the map to a database in the PC. The proces- 
sor employed cross-correlation to calculate the velocity 
vectors. A total of 800 sets of data was taken at each lo- 
cation for a specified Reynolds number (Re; i.e., the ratio 
of inertial forces to viscous forces). The selection of 800 
data sets was based on the examination of the data con- 
vergence. One set of data consisted of five PIV vector 
data for a 32 x 32 pixel interrogation area. These data 
were statistically averaged, and the mean vector fields 
were obtained and used for the examination of the flow 
structure. The measurements were performed in a clean 
room at the University Microsystem Laboratory at a 
controlled ambient temperature of 298 K. 

Methodology used (for electrophoretic mobility of DNA 
molecules and buffer solution EOF velocity) and 
temperature visualization 

Following [7], the electrophoretic velocity of the stained 
DNA molecules in the untreated PDMS channel with 
negligible electro-osmotic mobility was measured using 
(iPIV measurements. The total velocity of the seed parti- 
cles (i.e., DNA molecules) can also be measured through 
the (iPIV measurements for treated PDMS channels. 
With these velocities found, the bulk averaged EOF vel- 
ocity of the fluid {u) could be obtained following equa- 
tion (1) below: 

U= Um- Uph (1) 

where is the total velocity of the seed particles (i.e., 
DNA molecules) by (iPIV in treated PDMS channels, 
and Uph is the electrophoretic velocity of the DNA mole- 
cules in the untreated PDMS channel. 

With respect to measurement uncertainties, the most 
significant source of error was expected to be the mea- 
surements at the wall, and the biggest physical error in 
the (iPIV data was the Brownian diffusion of the stained 
DNA molecules. Out-of-plane Brownian diffusion causes 
a reduction of the signal-to-noise ratio of the cross- 
correlation peak, and such an error was estimated. Er- 
rors due to in-plane Brownian diffusion were essentially 
eliminated by temporally averaging the results in the 
steady flow. In fact, experimental errors due to the limit- 
ing spatial resolution of the CCD camera, as well as er- 
rors in determining magnification, were therefore the 
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major source of error in these results and found to be 
within ±15%. 

Visualization of the local fluid temperature was achieved 
with the same apparatus used for flow visualization and 
measurements (see Figure 3). Instead of using fluorescent 
particles, however, the channel was filled with a solution 
of rhodamine B, a fluorescent dye which shows a 
temperature-sensitive quantum yield in the range of 0°C 
to 100°C [5]. Experiments were conducted with a fluores- 
cence microscope equipped with a long-working dis- 
tance xlO objective lens. The images were recorded 
with the same equipment used for the (iPIV measure- 
ments. From the captured images, the detailed temperature 
distribution could be extracted. The intensity values of the 
captured images were converted to temperature using the 
intensity-versus-temperature calibration [5]. A calibration of 
the intensity of temperature was made for each solution. 

Sample preparation 

On the basis of standard lithography techniques, we 
constructed a 30-mm-long, 400-(im-wide, and 100- (im 
-high PDMS microchannel with a sudden contraction/ 
expansion (a ratio of 8:1:8) test section 20 mm in length. 
Reservoirs (4x4 mm) were cut at each end of the 
curved PDMS microchannel with a scalpel, and the 
channels were soaked for 12 h at 45°C in Ix TBE (Ix 
TEE contains, in 1 1, 108 g of Tris base, 55 g of boric 
acid, and 40 ml of 0.5 M EDTA, pH 8.3) to eliminate 
permeation-driven flow [3]. 

X-phage double -strand DNA (dsDNA) from New Eng- 
land Biolabs (Ipswich, MA, USA) was used as the tracer 
in the present study. The DNA was stained, with respect 
to the backbone, with a fluorescent dye (YOYO-1, 4.7:1 
bp/dye molecule), for a total length of 48.5 kbp DNA 
molecules, and diluted in Ix TBE. The dyed X-DNA had 
a contour length {1^) of 21 [im [3], and the longest relax- 
ation time (je) of 0.6 s (from uncoiled maximum length 
to coiled state) was measured and found in the present 
study. 

Results and discussion 

DNA molecule velocity profile with/without temperature 
effect 

Spanwise velocity profiles of DNA molecules at j = 0 in 
Ix TBE buffer at the inlet regions {x = 14.5 mm) of the 
Dh = 160 [im microchannel at = 5, 7.5, and 10 kV/m 
without joule heating are given in Figure 4a. The plug- 
like motion, a characteristic of an electrokinetic-driven 
flow, was apparent, and the velocity profiles remained 
fairly flat right to the wall for < 10 kV/m. On the 
other hand, the streamwise velocity profiles (not shown) 
of DNA molecules along the downstream at the inlet re- 
gime of the channel exhibited a nearly mountain-like 
distribution, similar to those reported in [3] for EOF 
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Figure 4 DNA molecule velocity at different heating 
temperatures and electric strength at the channel inlet, (a) 

Spanwise (x = 14.5 mm) and (b) transverse [x = 14.5 mm). 



with different magnitudes. The differences of about one 
order of magnitude were due to the former being elec- 
trokinetic driven, while the latter was pressure driven. In 
addition, the former was for DNA molecules along the 
downstream velocity, while the latter was for the EOF 
velocity of the buffer solutions. Nonetheless, they had 
the same developing trend, and they all increased as the 
£x increased. Figure 4b shows the corresponding trans- 
verse velocity distribution. Likewise, the similar plug/ 
uniform velocity profile again appeared. The insets in 
Figure 4a,b were made for clarity. Although the plug/ 
uniform velocity distribution in the y and z directions 
was what one would expect without the joule heating ef- 
fect, very small velocity differences in both the y and z 
directions were still noted upon close examination as 
the buffer solution was heated to different temperatures 
of 25°C, 35°C, 45°C, and 55°C. In addition, the velocity 
discrepancy increased as the heating temperature in- 
creased in both the y and z directions. 

Taking a further look at Figure 4a,b, with a working 
temperature of 25°C, the DNA molecule velocity had an 
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electrophoretic velocity apparently of the same order of 
magnitude in both the y and z directions due to the uni- 
formity across the stream. The velocity was constant 
across the channel width (along the z direction) and 
height (along the y direction) with an increase from 60 
to 125 [im/s as £x = 5 kV/m increased to 10 kV/m, as 
shown in Figure 4a,b. Tables 2 and 3 list the local velocity 
distributions measured at different electric strengths and 
heating temperatures, respectively. 

Figure 5 shows the relative velocity (|AV|; absolute 
value was taken), using a treated PDMS device to get the 
velocity of EOF of the buffer solution convection ob- 
served at four different temperatures, 25°C, 35°C, 45°C, 
and 55°C, for four corresponding heating powers at 
three electric strengths. Convection rates were estimated 
by the increases in buffer solution temperature. Two dif- 
ferent trends were observed: one (left half) at the same 
inlet position with different elevation and the other 
(right half) at the same elevation with a different down- 
stream position. The former showed an irregular velocity 
|AV| distribution as the heating temperature increased at 
different electric strengths, while the latter exhibited a def- 
inite quadratic |AV| increase as the heating temperature 
increased. The significant influence of the buffer solution 
temperature increase on DNA molecule stretching was 
clearly noted. This was partly because the rise in 
temperature of the buffer governed the evaporation rate 
which would result in suspension deposited on the chan- 
nel bottom surface, and consequently, it created the con- 
vective flux that dragged the DNA molecules and resulted 
in molecule stretching. Therefore, we concluded that ther- 
mal convection or thermophoresis caused DNA molecule 
elongation. In addition, with such temperature differences 
in the buffer solutions, certain effects due to thermal diffu- 
sion might have also occurred, which will be discussed 
later. 

Again, Figure 5 shows the velocity of the buffer solution 
convection observed for four different heating tempera- 
tures at the up, middle, and downstream locations, re- 
spectively (right half). The convection rates were 
approximately linear with the heating power and coin- 
cided with those found in Mao et al. [8], but they were 
strongly affected by the location where the velocity was 
measured. It was found that the convection effect became 



more dominant as the flow proceeded downstream, which 
was in good agreement with those of the temperature dis- 
tributions, namely, the temperature gradient became 
steeper downstream than upstream. 

DNA electrophoretic mobility and diffusion coefficient 

Electrophoresis is the net migration of a molecule induced 
by Coulomb forces on a charged molecule or particle. 
Despite the complexity of the physics that governs DNA 
electrophoresis, based on the above-stated velocity results, 
the electrophoretic mobility of long DNA in the buffers 
was found to be in the range of /u^p = 1.25 x 10"^ m^/Vs, 
which was in good agreement at a same order (approxi- 
mately 10~^) with [9]. Note that the thermophoresis effect 
in the calculation was neglected here for simplicity. 
Figure 6a shows the electrophoretic mobility of the DNA 
molecules. Generally, distribution is a linear function of a 
velocity-versus-electric field strength graph. In this figure, 
the slope of the lines represents the electrophoretic mobil- 
ity, with a close-up view of ^ at different temperatures. 
The temperature effect is not clearly noted. Again, this in- 
dicates that thermophoresis can be neglected. Further- 
more, the results from [10] were with ssDNA, which has a 
smaller molecular weight than the DNA molecules used 
in the present study. Thus, there was a much higher mo- 
bility of //ph , as depicted in Figure 6a. 

Diffusion in the present study could be classified as 
translational diffusion or rotational diffusion. Only trans- 
lational diffusion, i.e., diffusion of the center of the mass 
of DNA molecules, was considered. The translational 
diffusion was proportional to the thermal energy and, 
thus, proportional to k^T, as weU as the effective viscous 
mobflity, /u, FoUowing [8], we approximated the DNA 
diffusion coefficient as D = k^T I6jirjr^, where rj is the 
buffer viscosity, Tg is the gyration radius of the DNA 
molecule, and k^ is the Boltzmann constant, which was 
3.12 to 3.43 X 10"^ (im^/s in the temperature range of 
25°C to 55°C, as shown in Figure 6b. Further compari- 
sons were made with those of previous studies for //gp 
and diffusion coefficient A and the results are shown in 
Figure 6a,b, respectively. Given the different buffer solu- 
tions at different temperatures and the shorter gyration 
radius of the present study, as expected, the diffusion co- 
efficient D was lower, as iUustrated in Figure 6b. 



Table 2 Local velocity map (|Jim/s) with different heating temperatures 



y (Mm) 




25°C 






35°C 






45°C 






55°C 




10 


62.51 


93.40 


124.45 


61.55 


94.23 


123.59 


62.55 


95.88 


124.79 


63.89 


95.46 


126.97 


0 


62.23 


93.33 


124.09 


61.83 


93.53 


123.57 


62.22 


94.29 


125.06 


63.74 


94.89 


126.57 


-10 


62.10 


93.30 


123.88 


62.61 


94.59 


124.68 


63.48 


93.98 


125.68 


63.35 


94.79 


126.30 


Error (%) 


0.66 


0.11 


0.46 


1.72 


1.13 


0.9 


2.03 


2.02 


0.71 


0.85 


0.71 


0.53 



Velocity map at different heights of the channel in x= 14.5 mm and y= 10 to -10 |jm. 
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Table 3 Local velocity map (|am/s) with different heating temperatures and electric fields 
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Velocity map at different locations of the channel in x= 14.5, 14.6, and 14.7 mm and y-0. 



Heating effect on DNA molecule stretching 

Using detailed (iLIF observations, thermophoresis, often 
called the Ludwig-Soret effect (thermal diffusion), was 
considered [14]. The investigation of the Soret effect in 
the buffer solution was based on the determination of 



the following transport coefficient: Dmd? mutual diffusion 
coefficient; D^, thermal diffusion coefficient; and 5't, 
Soret coefficient. Detailed calculation of the values of 
the above-stated parameters improved our basic under- 
standing of the exact stretching mechanisms involved in 
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Figure 5 Relative velocity of the buffer solution convection. Velocity gradient at different electric fields and at a definite channel inlet x = 
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Figure 6 DNA molecule mobility and diffusion coefficient 
distribution, (a) DNA electrophoresis velocity versus electric field 
and (b) relationship of diffusion coefficient and buffer solution 
temperatures [11-13]. 



this study. However, due to the limitation of the mea- 
surements, several physical quantities above were not 
available at this stage. Further study could include this 
aspect. Nevertheless, thermal convection, as well as dif- 
fusion, was still noted. 

Figure 7 shows these results at different streamwise 
electrical strengths without the joule effect (<10 kV/m) 
at different temperatures. Note that thermal expansion 
occurred at = 0. There were two groups with a similar 
developing tendency but different rates of increase: one 
at a heated temperature between 25°C and 35°C and the 
other between 35°C and 55°C, with two different slopes. 
Obviously, the latter had a greater heating effect than 
the former as far as the stretching length was concerned. 
For all the electric strengths studied, the trend of the de- 
velopment of stretching versus temperature appeared to 
be similar. After deducting the thermal expansion length, 
the DNA molecule average stretching lengths were found, 
and they were plotted against applied electric fields, as 
shown in Figure 8. The most significant stretching hap- 
pened at £x = 10 kV/m as the heating temperature in- 
creased from 35°C to 55°C. The effect of electric strength 
that deducted the thermal effects was also as expected, al- 
though the rate of increase was minimal. As stated previ- 
ously. Figure 8 also shows the thermal expansion 
distribution {E^ = 0 kV/m) with different buffer tempera- 
tures. In addition, it was apparent that after the 
temperature rose to 45°C, the DNA molecule thermal ex- 
pansion coefficients appeared to be independent of 
temperature and reached a constant at about 0.097 K"^. 

Heating effect on the histogram of DNA stretch ratio 

Figure 9 shows the DNA histogram of the stretch ratio 
without the electric field applied at the inlet region. The 
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14.9 mm) via CLSM. 
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Figure 8 Average stretching length). After deducting the thermal expansion effect and coefficient of DNA thermal expansion versus 
temperature at the inlet region (14.6 to 14.9 mm) with the scale bar of 10 |jm.; 



heating effect was clearly noted as the maximum exten- 
sion length went from about 2.5 [im at 25°C to 6.5 [im 
at 55°C. In addition, 85% of the DNA molecules (^85%) 
were at 1.5 [im at 25°C versus 40% at 5.5 (im, even with 
no external electric field employed. The stretching was 
partly due to thermal expansion of the DNA molecules 
(<10%) and partly because of thermophoresis (>90%). 
Each contribution (10% versus 90%) can be calculated 
based on a measured thermal expansion coefficient in 
Figure 8 and obtained. 



Moreover, when electric strength was applied, the stretch 
ratio was enhanced. Figure 10 shows respectively the corre- 
sponding results at different regions (inlet/middle) with dif- 
ferent temperatures at = 10 kV/m and Deborah number 
(De) = 2.3. The effect of the position either at the inlet/or 
middle region can be seen. At the downstream middle re- 
gion, the DNA molecules seemed to be further stretched, 
and most significantly, more DNA molecules were found 
at a larger stretch ratio, for instance, 10% (inlet) versus 20% 
(middle) at 55°C and De = 2.3 for a stretch ratio of 0.4. 
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Figure 10 Histogram of thie stretcli ratio of DNA molecule after deducting the thermal expansion effect. At Ey^= ]0 kV/m at different 
temperatures. Inlet region: (a) 25°C, (b) 35°C, (c) 45°C, and (d) 55°C. Middle region: (e) 25°C, (f) 35°C, (g) 45°C, and (U) 55°C 



Stretching force distribution 

Extracting the data from Figure 10, the maximum exten- 
sion distribution was deduced to be a function of the 
stretching force. The stretching portions of the force- 
extension curves as a function of temperature are shown 
in Figure 11, in which the DNA molecule maximum ex- 
tension length versus hydrodynamic force after deducting 
the thermal effect can be drawn and compared with those 
from the well-known force law of the wormlike-chain 
(WLC) model. The stretching force clearly decreased as 
the temperature increased due to thermal convection and/ 
or thermophoresis, as evidenced by the thermal convec- 
tion velocity distributions, as shown in Figure 4b and 



especially in Figure 5a,b,c,d,e,f. With the thermal expan- 
sion effect deducted, the different temperature results 
were shown in Figure 11a. As expected, the temperature 
effect had a significant influence on extension. Unlike 
those in Hsieh et al. [2] or Hsieh and Liou [3], the present 
stretching behavior at a temperature of 55°C changed fol- 
lowing the evolution of double strand, transition, and sin- 
gle strand, based on CLSM in situ observation. Even so, 
similar linear dependence behavior was still found with 
different slopes. In fact, as the temperature increased from 
25°C to 45°C, the force-stretching length relationship was 
similar to that in Hsieh and Liou [3], as the comparison in 
Figure 11a shows. In addition, there was good agreement 
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Figure 1 1 Stretching portions of tlie force-extension curves as 


a function of temperature, (a) Maximum DNA molecule 




hydrodynamic force versus extension after deducting the thermal 




expansion effect, (b) Hydrodynamic force of the present study 




versus the force law from the WLC model. 





with almost equal slopes as the temperature increased 
from 25°C to 45°C. This finding also verified the correct- 
ness of the present measurement. 

As Figure lib shows, the present experimental data 
could be approximately fitted by applying the well- 
known WLC model. The hydrodynamic force measured 
and calculated through the Stokes formula was found to 
be a power law function of the forces via the WLC 
model with different exponents of 3.05 to 4.56 and dif- 
ferent coefficients (Ci to C4) with different temperatures. 
Obviously, the stretching forces were greater than those 
predicted by the WLC. Furthermore, the temperature ef- 
fect was again noted; otherwise the exponent found 
would have been nominally the same. 

Conclusions 

DNA molecule dynamics in gradual/sudden converging- 
diverging heated microchannels were extensively examined 
via CLSM visualization and (iPIV velocity measurements of 



dsDNA molecules in solutions at different tempera- 
tures, i.e., 25°C, 35°C, 45°C, and 55°C. The important 
points from this study were as follows: 

1. A decrease in the stretching force was observed as 
the solution temperature increased, which was in 
good agreement with that in Williams et al. [15]. 

2. Although thermophoretic stretching was not clearly 
noted, the effect still seemed to occur and to increase 
as the temperature increased. 

3. In addition to electrophoretic stretching, thermal 
convection made an equal contribution in terms of 
DNA molecule stretching. 

4. As a result of points 2 and 3, when the buffer 
solution temperature increased, the stretch ratio was 
three to four times that of the isothermal buffer 
solution. 

5. DNA molecule thermal expansion played a 
significant role (>50%) in DNA molecule stretching. 
Therefore, the present stretching mechanism 
included thermal expansion, thermal diffusion 
(thermophoresis), and thermal convection. 

6. Electrophoretic mobility and the translational 
diffusion coefficient of the DNA molecules were 
obtained and compared with those of existing data. 
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